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Q Context. Massive stars play a dominant role in the Universe, but one of the main drivers for their evolution, their mass loss, 
^ remains poorly understood. 

C/3 Aims. In this study, we test the theoretically predicted mass-loss behaviour as a function of stellar effective temperature across 
the so-called 'bi-stability' jump. 

^ , Methods. We observe OB supergiants in the spectral range 08-B3 at radio wavelengths to measure their thermal radio flux 
T""] ' densities, and complement these measurements with data from the literature. We derive the radio mass-loss rates and wind 
I^nJ ' efficiencies, and compare our results with Hq mass-loss rates and predictions based on radiation-driven wind models. 

Results. The wind efficiency shows the possible presence of a local maximum around an effective temperature of 21 000 K - in 
qualitative agreement with predictions. Furthermore, we find that the absolute values of the radio mass-loss rates show good 
agreement with empirical Ha rates derived assuming homogeneous winds - for the spectral range under consideration. However, 
the empirical mass-loss rates are larger (by a factor of a few) than the predicted rates from radiation-driven wind theory for 
objects above the bi-stability jump (BSJ) temperature, whilst they are smaller (by a factor of a few) for stars below the BSJ 
temperature. The reason for these discrepancies remains as yet unresolved. A new wind momenta-luminosity relation (WLR) 
for 08-BO stars has been derived using the radio observations. The validity of the WLR as a function of the fitting parameter 
related to the force multiplier QcH (Kudritzki & Puis 2000) is discussed. 

Conclusions. Our most interesting finding is that the qualitative behaviour of the empirical wind efficiencies with effective 
temperature is in line with the predicted behaviour, and this presents the first hint of empirical evidence for the predicted 
mass-loss bi-stability jump. However, a larger sample of stars around the BSJ needs to be observed to confirm this finding. 
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ABSTRACT 



1. Introduction 

Massive stars are the main drivers for the evolution of 
galaxies: they are the prime contributors to the energy and 
momentum input into the interstellar medium through 
stellar winds and supernovae, they generate the bulk of 
ionising radiation, and they are important sources for the 
chemical enrichment of carbon, oxygen and nitrogen in the 



Universe. Their very evolution through subsequent stages 
(OB -> LBV -> WR -> SN) is beheved to be driven by 
mass loss (e.g. Chiosi & Maeder 1986). Heger et al. (2003) 
have recently reviewed the way in which massive stars end 
their lives and they find that the type of remnant that is 
left - neutron star, black hole, or no remnant at all ~ is 
primarily dependent on the amount of mass lost via stellar 
winds. 
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The physical mechanism for these winds of massive 
OB stars has long been identified to be that of radiation 
pressure on millions of spectral lines (Lucy & Solomon 
1970; Castor Abbott fc Klein 1975: hereafter CAK). 
These radiation-driven wind models predict the mass-loss 
rates M of O stars to depend strongly on the stellar 
luminosity, and the terminal flow velocity (voo) to be 
proportional to the escape velocity (wosc), reaching wind 
velocities of thousands of km s~^. Over the last decades, 
the theory has proven to be very successful in explaining 
the overall mass-loss properties of O stars, with M of O 
stars up to 10~^ Mq yr~^ and a clear correlation between 
the wind terminal and the escape velocity: (woo/^^esc) — 
2.6 (Kudritzki & Puis 2000). 

Nevertheless, there are many open issues in the field 
of massive star evolution and their winds - even during 
the "best understood" O & B evolutionary phases. One 
of them is that of wind clumping (e.g. Eversberg et 
al. 1998). Many modern evolution models (Maedcr & 
Meynet 2003, Limongi & Ciofii 2006, Eldridge & Vink 
2006) use OB mass-loss rates from smooth wind Monte 
Carlo mass- loss predictions of Vink et al. (2000a, 2000b). 
However, the sizes of these mass-loss rates have recently 
been questioned by a number of ultraviolet (UV) studies, 
and downward revisions by factors of five to ten have 
been suggested (Bouret et al. 2003, FuUerton et al. 
2006). There are theoretical and observational indications 
(e.g. Evans et al. 2004) that the winds of OB stars are 
clumped, but what remains uncertain is by what factor 
the mass-loss rates could be over-predicted; is it as 
dramatic (i.e. factors ^ 10) as has been suggested by the 
UV studies, or is the clumping factor relatively small (a 
factor of two)?. It is relevant to note that a modest 
amount of wind clumping is actually required to match 
the Vink et al. (2000a, 2001) mass-loss predictions. If 
the clumping factor (/) is only 5 - the corresponding 
mass-loss rates arc a factor \/5 lower- (Repolust et al. 
2004, Mokiem et al. 2007), massive star evolution is not 
anticipated to be dramatically affected by clumping. 
For a recent discussion on wind clumping from optical 
studies, we refer to Puis et al. (2006), and Davies et al. 
(2007). 

Another topical issue is that of the "bi-stability jump" 
(BSJ): an empirical drop in the ratio (t^oo/^^esc) from 2.6 to 
1.3 around temperatures of about 21 000 K (Lamcrs et al. 
1995). Theoretically, the wind characteristics as a function 
of stellar spectral type are best described in terms of the 
wind efficiency number: r) = (Mvrx.)/{L^/c), a measure 
for how much momentum from the photons is transferred 
to the ions of the outflowing wind. Vink et al. (2000a) 
computed a large grid of wind models as a function of 
effective temperature and presented the following overall 
behaviour of the mass-loss rate as a function of decreasing 
temperature (Fig. 1). 

Over the temperature range between 50 000 K down 
to 27 500 K, the mass-loss rates drop rapidly. The reason 
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Fig. 1. The predicted wind efficiency number: rj = 
{Mvoo)l{L*/c) as a function of effective temperature 
(from Vink et al. 2000a). Note the presence of a local 
maximum at the position of the bi-stability jump around 
Teff = 21 OOOK. 



is that of a growing mismatch between the wavelengths 
of the maximum opacity in the UV and the flux maxi- 
mum, which progressively moves from extreme UV (Ly 
continuum) to far UV, from early O to early B stars. The 
behaviour is however reversed around the BSJ, where rj 
is now predicted to increase by a factor of 2-3 to a lo- 
cal maximum when Fe IV recombines to the more efficient 
Fc III configuration (Vink et al. 1999, Vink 2000). Below 
the BSJ, the first effect returns, and r] is again predicted 
to decrease with temperature. 

For B supergiants, Vink et al. (2000a) reported signifi- 
cant discrepancies between empirical and predicted mass- 
loss rates. Recent studies involving Ha (Trundle & Lennon 
2005, Crowther et al. 2006) and UV analyses (Prinja et al. 
2005) have reconfirmed these findings. However, discrep- 
ancies between empirical UV rates and radiation-driven 
wind models have now also been reported for O stars (sec 
FuUerton et al. 2006). It is currently unclear whether the 
reported discrepancies for B supergiants are due to the as- 
sumption of smooth winds (and the neglect of wind clump- 
ing) in the radiation-driven wind model for all spectral 
types, or instead that they are somehow related to the 
physics of the BSJ. 

1.1. The bi-stability jump 

The bi-stability jump was first discussed by Pauldrach & 
Puis (1990) in the context of their model calculations of 
an individual star, the Luminous Blue Variable (LBV) P 
Cygni. The first empirical work on the BSJ for a sample 
of stars was that by Lamers et al. (1995), who found a 
discontinuity in the ratio (uoo/i'esc) at spectral type Bl 
(Tofr ^ 21 000 K). For the earlier-type stars the ratio 
is ~ 2.6, but for the lower temperatures it drops to ~ 
1.3. They also discussed the possible existence of a sec- 
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ond BSJ at spectral type AO. We note that Crowther et 

al. (2006) performed a sophisticated analysis using line- 
blanketed model atmospheres (of Hillier & Miller 1998) 
of the temperatures of objects around the Bl BSJ and 
concluded the empirical BSJ in wind velocity represents a 
more gradual decrease than suggested by the step function 
of Lamers et al. (1995). 

On the theoretical front, Vink ct al. (1999) studied the 
nature of the winds around the BSJ and found a gradual 
change in wind properties, while they attributed the ori- 
gin of the BSJ to a gradual change in the ionization of the 
Fe lines that drive the wind. For a range of stellar param- 
eters and adopted velocity laws, Vink et al. (1999, 2000a) 
found an increase in the wind efficiency 77 by a factor 2-3 
independent of the assumed terminal velocity. When ac- 
counting for the empirical drop in terminal velocity, they 
found a jump in the mass-loss rate by a factor of five for 
stars at the same luminosity. They argued that the jump 
in mass loss is accompanied by a decrease of the ratio 
(■yoo/uesc), which is the observed bi-stability jump in ter- 
minal velocity. Using self-consistent models they obtained 
ratios within 10 percent of the observed values around the 
jump (Vink ct al. 1999, 2000b; Vink 2000) after a detailed 
investigation of the line acceleration for models around 
the jump, they demonstrated that M increases at the BSJ 
due to an increase in the line acceleration of Fe III near 
the stellar photosphere. 

The bi-stability mechanism may be the crucial piece of 
physics for understanding the least understood phases of 
massive star evolution, i.e. the LBV and the B[e] phases. 
Vink & de Koter (2002) found that the bi-stability jump 
may explain the mass- loss variability of LB Vs. Smith et 
al. (2004) showed that, under certain conditions (when 
stars have lost a lot of mass due to mass loss), the bi- 
stability jump may be able to explain the formation of 
pseudo-photospheres and possibly even the existence of 
Yellow Hypergiants. Kotak & Vink (2006) recently sug- 
gested that the quasi-sinusoidal modulations in the radio 
lightcurves of a number of supernovae may be caused by 
variable LBV wind strengths, suggesting that LBVs may 
be the progenitors of some supernovae, which would have 
profound consequences for our most basic understanding 
of massive star evolution. 

For rapidly rotating stars with radiative envelopes, the 
pole may be expected to be hotter than the stellar equator, 
as due to the Von Zeipel gravity darkening effect, imply- 
ing that the BSJ could provide a higher mass loss from 
the equator than from the pole for B stars. This will give 
rise to a high density "disk-forming region" at the stellar 
equator - possibly causing the B[e] phenomenon (Lamers 
& Pauldrach 1991, Pelupessy et al. 2000, Cure et al. 2005). 

However attractive the physics of the BSJ may be, 
the expected jump in mass loss has yet to be investigated 
on observational grounds, which is the prime goal of this 
study. We perform a comprehensive radio analysis of 
massive stars over the spectral range 08-B3, correspond- 
ing to Teff in the critical range of 35 000 - 15 000 K. 



The contents of the paper is as follows. Section 2 
briefly reviews the derivation of stellar mass-loss rates by 
means of radio data. Section 3 describes the stellar sam- 
ple, including the radio observations carried out specif- 
ically for this work. In Sect. 4 we explain the adopted 
stellar parameters and the derivation of radio mass-loss 
rates. Section 5 presents the wind efficiencies and a wind 
momentum-luminosity relation obtained with radio data. 
In Section 6, we compare our radio results with theoret- 
ical predictions. In Section 7 we analyze the correlation 
between radio- and Ha-derived mass-loss rates. Section 8 
closes with a summary and the conclusions of the present 
study. 

2. Stellar mass-loss rates from radio observations 

The measurement of stellar mass-loss rates can be 
achieved by means of at least three methods, which involve 
observations at radio, optical (Ha), and UV wavelengths. 
Each technique retrieves information from a different part 
of the wind. Although the radio method is not as sensitive 
as the Ha or UV methods, the radio flux density method 
is generally accepted to be the most accurate (Abbott et 
al. 1980, Lamers & Leitherer 1993). The reason is that the 
method does not require a priori knowledge of the excita- 
tion, ionization or velocity structure in the winds, as the 
other approaches do to a more or lesser extent. 

Although the radio method is considered to be the 
most accurate, it nonetheless depends on the assumption 
that the radio flux density is predominantly thermal in 
nature. This may be quite a good assumption for the ob- 
jects under consideration in this study, as we expect that 
non-thermal emission at these spectral types should not be 
important (even irrespective of the wavelength): all (13) 
targets in the range 08-B3 detected at two frequencies by 
Bieging et al. (1989), and by Scuderi et al. (1998) showed 
thermal spectral indices. Other potential inaccuracies in 
the radio mass-loss rates are due to uncertainties in the 
stellar distance d and only slightly on the value of the ter- 
minal velocity of the winds. We discuss the effect of wind 
clumping in the last Sections. 

Despite the weakness of the thermal radio emission 
from OB stellar winds (usually less of 1 mJy), the num- 
ber of closeby (< 3 kpc) targets is enough to perform a 
statistical study. The thermal radio flux density S*^ at a 
frequency from an optically thick stellar wind, is related 
to the stellar mass-loss rate M as (Wright & Barlow 1975, 
Panagia & FeUi 1975): 

M = 5.32 X 10-4 {S,y-' d-' r ^ ji ^ -1 _ 

Here /x is the mean molecular weight of the ions, Z the 

rms ionic charge, 7 the mean number of electrons per ion, 
and gi, the Gaunt factor, approximated by: 

= 9.77 [1 + 0.13 log{OAT^^\Zu)-% (2) 
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3. Observations at radio wavelengths 

In order to derive mass-loss rates from continuum radio 
observations, for stars with effective temperatures around 
21 000 K, we aimed at supergiants of spectral types 
in the range 08 to B3. We present here the results of 
our first observing campaign, with targets of spectral 
type spread over the mentioned range. The sample was 
completed with the addition of the results from previous 
observations taken from the literature and carried out by 
several authors. 

The observing campaign consisted of 12 targets with 
declinations between —57° to +37°. Eight of them were 
observed with the NRAC0 Very Large Array (VLA), 
while the southern ones were observed using the Australia 
Telescope Compact Array (ATCAj^. The stars can be 
identified in Table 1, with the acronyms ATCA or VLA 
(third column). The targets were required to be closeby 
and have as-accurate-as-possible determinations of termi- 
nal velocities. The integration time was estimated as the 
time needed to detect the radio flux density correspond- 
ing to the star mass-loss rate as predicted by Vink et al. 
recipell (Vink 2000). 

3.1. ATCA data 

The ATCA observations were conducted on June 26-27, 
2005, at 6B array configuration. The observing wave- 
lengths were chosen at the 12mm band, i.e. high frequen- 
cies, where not only the thermal flux density is stronger, 
but also the contribution of non-thermal emission to the 
radio flux density may be disregarded. The angular reso- 
lution is about 1". The observing frequencies were 17.728 
and 17.856 GHz, and the total bandwidth was 128 MHz 
in each of the two IFs. Due to bad weather 25% of the 
data were not useful, and had to be flagged. The sources 
were observed during 6 min scans, interleaved with 2 min 
scans of phase calibrator observations. The bandpass cal- 
ibrator was 1253-055, and the flux density calibrator was 
PKS 1934-638 (5'i2mm = 1-03 Jy) during both runs. Table 
1 lists the phase calibrator and the integration time for 
each target source. 

The data reduction and analysis were performed with 
the MIRIAD package. One star was detected, HD 148379, 
with a flux density of 0.28 ± 0.05 mJy (ure 2). The syn- 
thesized beam for each source is given in Table 1. Table 2 
shows the flux densities and l-cr errors, and, for the case of 
non-detections, the flux density value quoted as an upper 
limit corresponds to a 3-cr value. 



^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 

^ The Australia Telescope Compact Array is funded by 
the Commonwealth of Australia for operation as a National 
Facility by CSIRO. 

^ http://www.astro.keele.ac.uk/~jsv/ 



Table 1. Observations details 



Star 


Phase 


Tel. 


Int. 


Synth. 






Cal. 




time 
(h) 


beam 

(arcsec) 


HD 


42087 


0559-^238 


VLA 


0.81 


2.2 X 1.8 


HD 


43384 


0559-^238 


VLA 


0.83 


2.2 X 1.8 


HD 


47432 


0641-033 


VLA 


0.75 


2.5 X 2.3 


HD 


76341 


0828-375 


VLA 


0.50 


11.3 X 2.3 


HD 


112244 


j 1326-5256 


ATCA 


3.18 


0.8 X 0.4 


HD 


148379 


1600-44 


ATCA 


3.58 


1.0 X 0.8 


HD 


154090 


1626-298 


VLA 


0.75 


7.3 X 2.6 


HD 


156154 


1744-312 


VLA 


0.5 


7.8 X 2.6 


HD 


157246 


1657-56 


ATCA 


6.9 


0.6 X 0.5 


HD 


165024 


1740-517 


ATCA 


6.36 


0.7 X 0.5 


HD 


204172 


2052+365 


VLA 


1.25 


3.0 X 2.7 


BD- 


11 4586 


1832-105 


VLA 


0.5 


3.8 X 2.7 



3.2. VLA data 

The VLA targets were observed on September 9 and 17, 
2005, at C array, during a total observing time of 8 h. The 
flux density cahbrator was 1331+305 (S's.scm — 5.23 Jy). 
The integration time on source, synthesized beam, and 
phase calibrator corresponding to each source are given in 
Table 1. The observations were carried out at 3.5 cm with 
two IF pairs of 50 MHz bandwidth, to obtain a reasonable 
angular resolution with a moderate amount of calibration 
time, and lower noise receivers. 

The data was reduced and analyzed with the AiPS and 
MIRIAD routines. Two stars were detected above the noise: 
HD 76341, and HD 154090, with signal-to- noise ratio > 6 
(see Figure 2). The stellar field of HD 156154 showed weak 
emission at a 3-cr level, and a 2-a (probable) noise source 
coincident with the stellar position. Diffuse emission over 
the field of HD 42087 precluded the detection of the stellar 
wind, also superposed to a 2MASS source, and the infrared 
source IRAS 06066+2307. 

An extended source was discovered in the stellar 
field of BD-11 4586, centered at (i?A, Dec)j2ooo = 
(18''18'"6.744^-ll°18'16.83"), with a total flux density 
of 2.8 mJy. No correlation at any wavelength was found 
for this source in the literature. 

3.3. Previous radio observations 

We searched the literature for continuum radio observa- 
tions of late OB supergiants with spectral types between 
08 to B3. Among the contributions considered, the com- 
prehensive study by Bieging et al. (1989) on OB stars 
listed high-resolution radio observations of 16 such super- 
giants, six of which were detected. Lamers & Leitherer 
(1993) presented seven radio detections. Scuderi et al. 
(1998) obtained radio flux densities -or upper limits- for 
ten candidates. We have also included isolated results from 
Stevens (2006) (private communication), Benaglia et al. 
(2001), and the latest radio flux densities detected by Puis 
et al. (2006) on two 09.5 targets. From the older studies 
(Bieging et al. 1989, Lamers & Leitherer 1993), only detec- 
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Fig. 2. Left panel: VLA radio continuum image of HD 154090 at 8.46 GHz. Contour levels: -0.08, 0.08 (2cr), 0.12, 0.15, 
0.18. and 0.21 mJy beam^^. Central panel: VLA radio continuum image of HD 76341 at 8.46 GHz. Contour levels: 
-0.15, 0.15 {3ct), 0.2, 0.25, 0.3 and 0.35 mJy beam"!. Right panel: ATCA radio continuum image of HD 148379 at 
17.8 GHz. Contour levels: -0.1, 0.15 (3(7), 0.2, 0.25, and 0.3 mJy beam^^. The optical position of the star is marked 
with a cross, and the synthesized beam is displayed in the bottom left corner of each image. 



tions were included, as the observations were carried out 
using the former VLA less sensitive receivers. The star 
HD 193237 (observations by Scuderi et al. 1998) is to be 
discussed elsewhere amongst LBVs. The star HD 152408 
is not included within the present sample due to its am- 
biguous classification as either "08 lafpe" or "WN9 ha" 
(Bohannan & Crowther 1999). 

The star HD 149404 is listed by Lamers & Leitherer 
(1993). It has been studied by Rauw et al. (2001), and 
Thaller et al. (2001), through optical spectroscopy. They 
found evidences of a binary system, formed by an 07.5 1(f) 
and an ON9.7 I, with a period of ~ 9 d. Most of the radio 
emission would be produced in a colliding wind region; the 
results are consistent with a picture where the secondary 
star is undergoing Roche lobe overflow. Within this sce- 
nario, the radio flux density measured should come from 
the region of colliding winds, as well as the winds of both 
components. The mass-loss rate derived from the radio 
flux density is thus only an upper limit and we disregard 
from our upcoming analysis this object as well. 




32000 28000 24000 20000 
Teff (K) 



16000 



Fig. 3. Radio flux density corrected for the stellar dis- 
tance, as a function of effective temperature. The squares 
represent detected stars, whilst the triangles represent ra- 
dio flux density upper limits. The star HD 80077 was ex- 
cluded due to its large error in distance (^40%). 



3.4. Considerations on radio-detectability 

The 08 to B3 supergiant sample comprises 19 detections 
(up to 3 mJy) plus upper limits to 11 sources (typically 
below 0.3 mJy). With the aim of comparing the different 
radio flux densities measured towards the sample stars, we 
have extrapolated the values of Si, to a unique frequency 
of 14.95 GHz, assuming that the thermal flux behaves as 
S„ (X v°-'^ (Setia Gunawan et al. 2000, WiUiams 1996). In 
Figure 3 we have plotted the quantity Si, * (c?/lkpc)^ oc L 
as a function of the stellar effective temperature. We found 
no correlation of the corrected flux with stellar effective 
temperature, although the errors in stellar distances are 
strongly affecting the result. 

No correlations between radio flux density and termi- 
nal velocity, distance or visual magnitude were found, as 
expected. 



4. Stellar parameters and mass-loss rates 

Table 2 lists the 30 selected stars in two groups (observed 
by us, and compiled from the literature), ordered by HD 
number, together with other observational stellar param- 
eters. 

The spectral types for O stars were taken from the 
COS catalogue (Mai'z Apellaniz et al. 2004), Crowther et 
al. (2006), and Puis et al. (2006). The various references 
on spectral types of B stars can be seen in Table 2. The 
V, B and absolute adopted magnitudes are also given. 

Most of the stellar wind terminal velocities quoted 
were published by Howarth et al. (1997). For four objects, 
the values were taken from the interpolations by Prinja 
et al. (1990), their Table 3, and involved relatively large 
errors. For the rest, we have adopted errors of 10%. 

To achieve the goal of a uniform database, we have 
established a common temperature scale for 08-B3 su- 
pergiants, re-derived the stellar distances, and calculated 
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Table 2. Supergiants with radio observations: measure- 
ments 
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HD 151804 


08 laf^ 


LL93 


PS 22 5 26^ 


1445*^ 


-6.3-^ 


1.9" 


0.3 


4.800 


0.4 


0.1 


HD 152236 


B1.5 la^-'^ 


Ste06 


4.73, 5.21^* 


390^ 


-8.8^ 


2.0^ 


0.3 


8.640 


2.4 


0.1 


HD 152424 


OC9.7 la" 


LL93 


6.31, 6.71' 


1760^ 


-6.3^ 


1.7" 


0.4 


8.000 


0.21 


0.04 


HD 163181 


BN0.5 lap^ 


BCKOl 


6.60, 7.15"= 


520^ 


-6.1^ 


1.6" 


0.4 


8.640 


0.44 


0.05 


HD 169454 


Bl la^ 


BAC89 


6.65, 7.39"= 


850^ 


-6.8^ 


0.9" 


0.2 


15.00 


1.9 


0.1 


HD 190603 


B1.5 Ia+'^ 


SPS98 


5.65, 6.19'* 


485^ 


-7.5^ 


1.5^ 


0.4 


14.95 


< 0.7 




HD 194279 


B2 la'^ 


SPS98 


7.05, 8.07'* 


550^ 


-7.0^ 


1.2^ 


0.3 


8.450 


0.44 


0.1 


HD 195592 


09.5 


SPS98 


7.08, 7.95' 


1765'' 


-6.3^ 


1.0" 


0.2 


14.95 


0.9 


0.13 


HD 198478 


B2.5 la'^ 


SPS98 


4.86, 5.28'* 


470^ 


-6.4^ 


0.8^^ 


0.2 


8.450 


< 0.27 




HD 209975 


09.5 Ib^ 


PMS06 


5.10, 5.18' 


2050^ 


-5.45^ 


0.8^ 


0.2 


14.94 


n.422 


0.12 


CygOB2-10 


09.5 F 


PMS06 


9.88 11.47" 


1650^ 


-6.95^ 


1.7^ 


0.2 


14.94 


0.3 


0.1 


BD-14 5037 


B1.5 la^ 


SPS98 


8.22, 9.58"= 


750^^ 


-7.4° 


1.0" 


0.2 


14.95 


1.1 


0.2 



'^Perryman et al. 1997, 'Johnson et al. 1966, "Howarth et al. 1997, '*Crowther et al. 2006, "derived with CHORIZOS: see text, 
■^Lesh et al. 1968, "Moffett & Barnes 1979, ''Walborn 1976, 'GOS Catalogue (Mai'z Apellaniz et al. 2004), •'Martins et al. 
2005, '^Garrison et al. 1977, 'Walborn 1973, ™Hiltner et al. 1969, "Feinstein & Marraco 1979, "Prinja et al. 1990, "Crawford et 
al. 1971, "Walborn 1982, ''Walborn 1972a, ^Lcitlicrcr ct al. 1995 (LCK95), *Schild ct al. 1983, "Bohannan & Crowthcr 1999, 
"Humphreys 1978, ""Walborn 1972b, ^Kozok et al. 1985, ^Bieging et al. 1989 (BAC89), "^Puls et al. 2006 (PMS06), '"Massey & 
Johnson 1991, "Morgan et al. 1955. SPS98: Scuderi et al. 1998. LL93: Lamers & Leitherer 1993. Ste06: I.R. Stevens (private 
communication). BCKOl: Benaglia et al. 2001. 



the mass-loss rates for all stars in the sample using those 
consistent parameters. 

The effective temperature scale was established taking 
into account the latest work by Crowther et al. (2006) for 
individual cases whenever possible, and also from interpo- 
lation using their Table 4. For some O supergiants the val- 
ues were taken from the models by Martins et al. (2005). 
Isolated values arc from Lamers & Leitherer (1993) and 
Puis et al. (2006). The resulting scale is shown in Figure 4. 
We have assigned to the thirteen different spectral types 
and subtypes covered, from 08 to B3, a correlated number 
(from 2 to 14). Figure 4 plots the stellar effective temper- 
ature as a function of this integer number. Spectral types 



were used to derive temperatures for the stars in the sam- 
ple. The uncertainties introduced by the discrete nature 
of the spectral classification scheme, the current state of 
the art of hot-star atmosphere modeling, and the possible 
slight mis-classifications are expected to be of the order 
of 2000 K for late-0 stars and 1500 K for early-B stars. 
The effective temperature adopted for each star is listed 
in Table 3. 

References for the adopted stellar luminosities are 
quoted in Table 3. The effective temperature and luminos- 
ity of HD 42987, HD 148379, and HD 151804 have been 
derived specifically for this work by means of Ha profile 
fitting (see Section 7). 
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Fig. 4. The adopted stellar effective temperatures as a 
function of spectral type, for OB supergiants. The straight 
line represents the fit y = 37327(±387) - 1583(±44)a;, 
where x is an index running on the different spectral types 
(from 2 to 14), and y = T^s- 



Fig. 5. Mass loss rates derived from radio contin- 
uum observations, as a function of effective temperature 
(squares), and upper limits (open triangles). 

here. 



In order to derive the stellar masses, needed in turn to 
obtain predicted mass-loss rates (see Section 6), we used 
the calibration of surface gravity {g) with effective tem- 
perature published by Crowther et al. (2006) (see their 
Figure 1, which was built using SMC data presented by 
Trundle & Lennon 2004). The likely range in the stellar 
mass for each star can be derived allowing a deviation by 
up to 0.2 dex in \og{g). 

Distances were derived using CHORIZOS (Maiz 
Apellaniz 2004) by comparing the observed with 
BVJHKs photometry with low-gravity Kurucz atmo- 
spheres. The temperature for each star was taken from 
Table 3 and both £^4405-5495 and -R5495 (the monochro- 
matic equivalents to Eb-v and respectively, were left 
as free parameters. CHORizos is a Bayesian code that finds 
all the solutions for the general inverse photometric prob- 
lem in which N parameters (two in this case) are fitted 
to M independent colors (three in this case). After apply- 
ing CHORIZOS to obtain the unextinguished photometry, 
the spectroscopic distance was calculated using the val- 
ues of Mv in Table 2, assuming uncertainties between 0.5 
and 1.0 magnitudes depending on the spectral type. In all 
cases we found that the largest contribution to the final 
uncertainty in the distance, e^, is that of the uncertainty 
in Mv and not that of the photometry itself. 

We investigated the membership of the sample stars 
to stellar clusters or OB associations. For those stars 
belonging to stellar groups, we found that in most cases, 
our obtained individual distances agree well between 
errors, with those of the clusters/associations. For the 
star HD 151804, catalogued as a member of Sco OBI 
(Humphreys 1978, d = 1.91 kpc), the association distance 
was adopted. Very recently, the distance to HD 169454 
was established as 1.36 kpc (Hunter et al. 2006). For the 
sake of consistency, we have adopted the distance derived 



In order to derive the radio mass-loss rate values for 
previously observed stars, we have taken into account the 
radio flux density at the highest frequency whenever pos- 
sible (see Table 2). 

All mass-loss rates values were derived using Eq. (1), 
and assuming the measured flux density is thermal. This 
looks like a reasonable hypothesis, as all but one star (HD 
151804) were observed at > 8 GHz. Besides, the sam- 
ple stars that have been observed by Bieging et al. (1989) 
and Scuderi et al. (1998) at more than one frequency, have 
shown thermal emission. For all but one star we have 
adopted abundances of 1h = 0.85, Fho = 0.15, and as- 
sumed that at the temperatures considered, all H is ion- 
ized, while He is singly ionized. Consequently /i = 1.45, 
Z = 7 = 1. The differences in the mass-loss rate results 
in considering /i ~ 1.3 - 1.6 were lower than 12%. If the 
He is doubly ionized in O stars, the mass-loss rates de- 
rived here are underestimated by about 30%, although we 
find this unlikely given the stellar effective temperatures 
involved. For HD 151804, we adopted Y^e = 0.25. Errors 
in ji, Z, and 7 were taken as 10%; aj, = 0.5L*. Mass- 
loss rate values and the corresponding errors derived from 
error propagation in Eq. (1) are listed in Table 3, and 
plotted in Fig. 5. In the case of the non-detections, a 3-cr 
flux density upper limit was used, and the corresponding 
mass-loss rates are upper limits. 

5. Wind efficiencies 

In order to better compare radio results with results from 

other wavelengths or models, we proceeded to compute 
the wind eSiciencies corresponding to the radio mass-loss 
rates. The values are shown in Table 3, and in Figure 
6 with filled squares. At Tcs ~ 21.5 kK there are two 
(squares) stars with wind efficiencies well above those of 
their neighbours: HD 169454 and HD 2905. The star HD 
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Table 3. Adopted stellar parameters and results 



Star 


Sj). typo 


T„,i 


k)g(L«/L.O 


log{g) 




V 


M 












(K) 










( 


10 Al 


yr 


' ) 


HD 156154 


08 lab(f) 


33179" 


5.68" 


3.10 


19.2 




< 4.4 




3.2 




HD 151804 


08 laf 


29000'' 


5.88'' 


3.00 


41.4 


1.0 ±0.12 


10.4 


3.70 


3.2 


16.0** 


BD-U 4586 


08 Ib(f) 


33179" 


5.68" 


3.10 


19.2 




< 5.6 




3.2 


(2.6)= 


HD 112244 


08.5 lab(f) 


32274" 


5.65" 


3.10 


20.0 




< 3.0 




2.6 




HD 76341 


09 lb 


31368" 


5.61" 


3.05 


18.2 


1.3 ±0.17 


7.2 


2.51 


2.1 




CygOB2-10 


09.5 I 


29700^* 


5.82'* 


3.00 


32.8 


0.5 ±0.07 


4.4 


1.77 


2.7 




HD 195592 


09.5 I 


29000-'' 


5.60-' 


3.00 


21.7 


1.0 ±0.14 


4.8 


1.69 


1.2 


3.8'' 


HD 30614 


09.5 la 


29000-^ 


5.63-^ 


3.00 


23.3 


0.6 ± 0.08 


3.3 


1.22 


1.5 


5.0-'' 


HD 209975 


09.5 lb 


32000"* 


5.31'' 


3.10 


9.5 


1.2 ± 0.15 


2.4 


0.94 


0.7 


l.l"* 


HD 37742 


09.7 lb 


28500-'' 


5.57" 


2.95 


19.4 


0.4 ±0.05 


1.5 


0.49 


0.9 


(1.8)'' 


HD 47432 


09.7 lb 


28500-'^ 


5.57" 


2.96 


19.8 




< 3.1 




1.1 


(1.9)'' 


HD 152424 


OC9.7 la 


28500-^ 


5.57" 


2.96 


19.8 


1.0 ±0.14 


4.5 


1.65 


1.0 




HD 37128 


BO la 


27000-^ 


5.44-^ 


2.90 


15.9 


0.5 ± 0.07 


1.6 


0.54 


0.5 


2.3^ 


HD 204172 


BO lb 


27500-'' 


5.50-^ 


2.93 


18.2 




< 5.0 




0.7 




HD 163181 


BN0.5 lap 


26000-'' 


5.57" 


2.85 


22.4 


0.1 ±0.02 


2.1 


0.73 


2.5 




HD 154090 


BO. 7 la 


22500-^ 


5.48-^ 


2.60 


18.1 


0.2 ± 0.03 


1.3 


0.48 


3.2 


1.0-'^ 


HD 2905 


BC0.7 la 


21500-^ 


5.52-'' 


2.53 


20.3 


0.4 ±0.05 


2.5 


0.84 


2.3 


2.0-'^ 


HD 169454 


Bl la 


21500-'' 


5.45-'' 


2.53 


17.3 


0.5 ± 0.06 


3.4 


1.14 


2.9 




HD 157246 


Bl lb 


20800^ 


5.45-^' 


2.45 


16.4 




< 0.6 




3.4 




BD-14 5037 


B1.5 la 


20500-'^ 


5.60-'^ 


2.45 


24.5 


0.2 ± 0.04 


2.5 


1.18 


4.9 


(3.4)^' 


HD 152236 


B1.5 Ia+ 


18000-^ 


6.10-^ 


2.20 


73.4 


0.1 ±0.02 


8.3 


2.70 


41.7 


6^ 


HD 190603 


B1.5 Ia+ 


18500-^ 


5.57-'' 


2.25 


21.8 




< 2.1 




6.8 


2.5^ 


HD 148379 


B1.5 Jape 


18500'' 


5.56'' 


2.25 


21.3 


0.1 ±0.01 


0.8 


0.29 


5.8 


1.5** 


HD 194279 


B2 la 


19000-^ 


5.37-^ 


2.30 


13.9 


0.2 ± 0.02 


1.5 


0.58 


3.3 


1.1^ 


HD 41117 


B2 la 


19000-^ 


5.65-'' 


2.30 


26.4 


0.1 ±0.02 


2.1 


0.76 


8.5 


1.0-'^ 


HD 80077 


B2 Ia+ 


18500-'' 


5.40-^' 


2.25 


14.7 


0.1 ±0.01 


1.7 


0.69 


30.2 




HD 165024 


B2 lb 


18000-'' 


5.40-' 


2.20 


14.7 




< 0.6 




0.9 




HD 198478 


B2.5 la 


16500-'^ 


5.03-'' 


2.05 


6.3 




< 0.5 




1.1 


0.2-*^ 


HD 42087 


B2.5 lb 


19000'' 


5.08'' 


2.30 


7.1 




< 1.1 




0.9 


0.2'' 


HD 43384 


B3 lab 


15500-^ 


4.88'' 


1.95 


4.4 




< 1.7 




0.3 


(0.3)'' 



"Martins et al. 2005, ''this work, '^Scuderi et al. 1992, ''Puis et al. 2006, ^Herrero et al. 2002, -''Crowther et al. 2006, ''Lamers & 
Leitherer 1993, ''Morel et al. 2004, 'Lamers et al. 1995, -'Scuderi et al. 1998. 



169454 has been represented twice, for wind efficiencies 
derived according to two possible distances (sec Sect. 4): 
with a square when the distance adopted here is used, and 
with a triangle when the Hunter et al. (2006) value is used. 

We wish to check if the empirical shapes of the wind 
efficiency could be an artifact of the fact that we observe 
objects with different masses and luminosities. This nec- 
essarily introduces a scatter in 77, but we may be able to 
correct for this. The correction was derived in the follow- 
ing manner. From models for different parameter sets of 
(M*, I/*), we obtain the corresponding wind efficiencies. 
We choose a parameter set for which 77 = 1 for the max- 
imum effective temperature and Woc/'i'csc ratio (sec Fig. 6 
from Vink et al. 2000a). The ij values were plotted as a 
function of their Eddington factor (cx M*/L*), resulting 
in Eddington-corrected rj values (Fig. 7). For each tar- 
get star we can compute the r] correction value according 
to its Eddington factor, and divide the uncorrected 77 by 
the correction value. The results are also shown in Figure 
6. Though an additional source of scatter is being intro- 
duced as a result of this via the errors in stellar masses 



(not shown in the Figure), the 'corrected' wind efficiencies 
are similar to the original values. Consequently, hereafter 
we will use non-corrected values. 

5.1. The empirical wind-momentum relation 

The radio observations allowed us to estimate the mass- 
loss rates for the 19 detected stars in a uniform way. 
We have utilised this database to provide a radio de- 
rived calibration for the coefficients Dq.x of the wind 
momenta-luminosity relation (WLR; see Kudritzki & Puis 
2000) for late O (08-09) and BO supergiants. Our WLR 
is shown in Figure 8. If the wind momenta equal D = 
M* (i?/R0)i/2, and log(L') = L>o + a; log(L/L0), a lin- 
ear fit yields the values Dq = 26.6 ± 2.1; a; = 1.0 ± 0.4 
when only the O stars (08-09) arc included, and Dq = 
25.0 ± 2.7; x = 1.3 ± 0.5 when the BO stars are also in- 
cluded in the fit. For obvious reasons Bl stars are excluded 
as they are located at the BSJ, whilst for spectral types 
later than Bl, we have too few datapoints to derive a re- 
fiable WLR. 
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Fig. 6. Observed wind efficiency versus effective temper- 
ature (squares). Crosses represent 'corrected' values (see 
text). The triangle corresponds to the star HD 169454 if 
the Hunter ct al. (2006) distance is considered. Note the 
possible presence of a local maximum at the predicted po- 
sition of the BSJ. 
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Fig. 7. Wind efficiencies taken for the twelve models from 
Fig. 5 from Vink et al. (2000a), as a function of the 
Eddington factor F. The straight line represents unity. 

We note the shallow slope of our derived WLR. For 
early-type O stars the slope, x, is significantly larger. In 
standard radiation-driven wind theory, the slope, x, is the 
reciprocal value of the "effective" CAK force multiplier 
aeff (see Kudritzki & Puis for definitions). For early O- 
stars, this value is approximately 2/3, the reason for which 
the stellar mass cancels and the WLR can be used in the 
first place. 

Here, we find an empirical aes in the range 0.78-1 for 

08-BO stars, which may be consistent with the high a 
value of a^^'" = 0.85 predicted through Monte Carlo simu- 
lations by Vink et al. (1999) for models at T^s = 30 OOOK. 
A high value of may suggest a relatively large value for 
the ratio of {voo/vesc)- Interestingly objects in this tem- 
perature range indeed show the largest (Voc/veac) ratio ac- 
cording to Fig. 6 of Kudritzki & Puis (2000) . We note how- 
ever that for these large values of a^s (when a^s no longer 
equals 2/3), the concept of the WLR may lose its meaning. 



Fig. 8. Wind momenta D of galactic late 0-, and early 
B supergiants as a function of luminosity, derived using 
radio continuum observations. 

Investigations of the wind efficiency r] versus temperature 
may therefore be a more valuable tool to understand wind 
driving for different effective temperature regimes. 

6. Mass-loss rates: expected versus predicted 

Once we have estimated the wind efficiencies from mass- 
loss rates derived using radio flux densities, we proceeded 
to compare the results with model predictions. Vink 
(2000), and Vink ct al. (2GGGa) developed mass- loss rate 
models and presented them for a (I/*,M*) grid. Based 
on the models, they produced a recipe to derive a stel- 
lar mass-loss rate if the stellar mass, luminosity, effective 
temperature and terminal velocity are known. We have 
made use of the routine to derive the expected mass-loss 
rates for the radio observed stars. The program also pro- 
vides information on the position of the star with respect 
to the BSJ, according to the specific stellar parameters. 
The predicted mass-loss rates Mprcd are given in Table 3. 

The recipe to derive the predicted mass-loss rates has 
helped to discriminate which stars are above and below the 
BSJ temperature, and we divided them into two groups. 
Both groups consist of 15 stars. The first group contains: 
HD 30614, HD 37128, HD 37742, HD 47432, HD 76341, 
HD 112244, HD 152424, HD 156154, HD 151804, HD 
163181, HD 195592, HD 209975, HD 204172, BD-11 4586, 
and Cyg 0B2 No. 10. The second group contains the re- 
maining 15 objects. 

The errors involved in the derivation of the predicted 
mass-loss rates are aL, = 0.5 i*, 0.2 dex in log(5'), and 
(TToff = 2000 — 1500 K. A direct comparison between ex- 
pected and measured values is presented in Figure 9. We 
show the wind efficiencies derived using radio observations 
towards Galactic stars, and, superposed, the predicted 
wind efficiency number constructed for log(L*/L0) = 6.0, 

= 60 Mq, Woo/'f^esc = 2.6, as in Fig. 1. 

Our results show that the stars above the BSJ 
temperature have predicted mass-loss rates lower than 
the empirical radio mass-loss rates. The reverse holds 
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Fig. 9. Modelled wind efficiency for log(L/L0) = 6.0, 
Af, — 60 M0, Vao/vcsc — 2.6 (dashed line, same as in Fig. 
1) superposed to the wind efficiencies derived by means of 
radio observations, for Galactic stars (squares). 



true for the stars below the BSJ temperature. If the radio 
mass-loss rates are affected by significant wind clumping, 
these empirical rates may need to be reduced by a certain 
amount, and possibly be brought into agreement with the 
predictions of radiation-driven wind models. If for stars 
above the BSJ temperature, the discrepancy between the 
empirical rates and the theoretical rates is attributable 
to wind clumping overestimating empirical rates - whilst 
the theoretical rates are not affected - then, adopting 
a similar clumping correction for stars below the BSJ 
temperature, would worsen the discrepancy between the 
empirical and theoretical rates for B supergiants. 

The reader should be aware that our analysis has been 
performed using derived mass-loss rates only. "Upper lim- 
its" are not taken into account in the comparison between 
empirical radio rates and mass- loss predictions. Therefore 
our comparison on the absolute value of the mass-loss rate 
is necessarily incomplete. Having noted this, there is no 
a priori reason why incompleteness should be a function 
of stellar temperature, and although the issue of the ab- 
solute mass-loss rates of hot-star winds is still an open 
one, the qualitative comparison between the empirical ra- 
dio mass-loss rates and the Monte Carlo predictions of 
radiation-driven wind models appears to support the pre- 
dicted mass-loss BSJ. More data around the BSJ will be 
necessary to confirm this. 



7. Comparison of radio and l-la results 

Many of the stars composing our sample have mass-loss 
rates in common with measurements from Ha profiles. The 
values of Mho are listed in Table 3, with their references 
(Lamers & Leitherer 1993; Scuderi et al. 1992, 1998; Puis 
et al. 2006; Crowther et al. 2006; Morel et al. 2004; Herrero 
et al. 2002). 




Fig. 10. Ha proffie of HD 148379. Source: UVES PoP 
Survey. 
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Fig. 11. Wind efficiencies from Ha measurements for 
the sample stars. Filled squares: results from CMF- 
GEn/fastwind profile fitting. Open squares: results from 
other methods. Dashed line: same as Figure 1. 



Additionally, three new Ha mass-loss rates have been 
derived specifically for this work. The stars HD 42087, 
HD 148379, and HD 151804 were analyzed in an identical 
manner to Crowther et al. (2006). For the first two stars 
the analysis was based on VLT/UVES data (PoP Survey, 
Bagnulo et al. 2003). For HD 151804, AAT/UCLES spec- 
troscopy was used, and the data was re-analyzed using 
current version of code plus TLUSTY velocity structure 
at depth. Figure 10 displays the Ha profile of HD 148379 
with errors, showing the reliability with which mass-loss 
rates can be achieved. 

Figure 11 presents the wind efficiencies derived us- 
ing mass-loss rates from Ha profiles for 20 stars out of 
the present sample (see Table 3). Results from for mod- 
ern CMFGEN and FASTWIND profile fitting are represented 
with filled squares. The rest (open squares) are specified 
in Table 3 with brackets. The results with Tcs are more 
scattered than they were for the radio wind efficiencies, 
but they do not contradict the predicted behaviour. 
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Fig. 12. Wind efRciencies from Ha measurements for 
SMC stars (from Trundle et al. 2004, and Trundle & 
Lennon 2005). 



An analogous plot to the Galactic Ha one can be con- 
structed for extragalactic stars. Trundle et al. (2004), and 
Trundle & Lennon (2005) have obtained the mass-loss 
rates of 16 B supergiants earlier than B4 in the Small 
Magellanic Cloud, by means of Ha measurements, all 
treated in a homogeneous way. These Mh„ are shown in 
Figure 12, as a function of the temperature scale adopted 
in this work. Possibly, the wind efficiencies for some stars 
around 22 000 K are above the overall decreasing trend, 
but larger samples are required to test if the SMC mass- 
loss rates indeed have a local maximum around the BSJ. 

Studies of the comparison between mass-loss rates 
derived using different methods have been conducted 
by various groups. Very recently, FuUerton et al. (2006) 
derived mass-loss rates from UV profiles and investigated 
whether the UV rates matched those estimated from 
Ha profiles and radio observations. They found that 
the results follow the relation: Afuv < Afna < -^j-adio 
and attributed this discrepancy to the presence of wind 
clumping. 

In the present study we have assumed that the free- free 
radio emission is due to smooth winds. If the winds are 
significantly clumped, the radio-rates derived in our pa- 
per are likely overestimates (see e.g. Blomme et al. 2003). 
We constructed Figure 13 with all radio and Ha mass- 
loss rate results available. From the figure it is clear that 
there is a good match between both sets of data, although 
the Ha lines are formed in the inner wind, whilst the ra- 
dio emission arises from the outer wind. This agreement 
does not necessarily discard clumping, but it may support 
the hypothesis that the clumping structure is uniformly 
distributed throughout the whole extension of the wind. 

A similar agreement between the inner and the outer 
wind was found for a large sample of O stars by Lamers & 
Leitherer (1993). However, these results have been chal- 
lenged by Fullerton et al. (2006) and Puis et al. (2006). 




log(M-Haipha) 



Fig. 13. Ha-derived versus radio-derived mass-loss rates. 
Units are 10~^ M0 yr~^. Squares: stars above the BSJ 
temperature; open circles: stars below the jump tem- 
perature. Symbols with arrows represent upper limits in 

A^fradio- 



8. Summary and conclusions 

Throughout the present investigation, a database of radio 
data of 30 08-B3 supergiants has been constructed, with 
observations carried out by different authors. The stellar 
distances have been reviewed and re-derived, using a sub- 
routine that takes into account the IR and optical indices. 
We have defined a uniform stellar temperature scale for 
the spectral range 08-B3. 

We have performed high resolution continuum ra- 
dio observations towards 12 08-B3 supergiants, detecting 
three of them: HD 76341, HD 148379 and HD 154090, and 
derived their mass-loss rates. The rates of another 18 stars 
were re-determined using results from the literature. 

New Ha mass-loss rates have been obtained for three 
objects: HD 42087, HD 148379, and HD 151804. The 
comparison of mass-loss rates derived from Ha and ra- 
dio observations yielded rather good agreement, contrary 
to some results from previous investigations. However, we 
note that there are outliers/significant discrepancies for 
individual cases. 

The renormalization of radio observations allowed us 
to derive a new calibration for the WLR for supergiants 
of spectral types in the range 08 to BO. 

The comparison between predicted (p-values) and 
radio-derived (r-values) mass-loss rates yields different 
results for stars above and below the BSJ temperature. 
For the latter group we find p-values that are larger 
than the r-values, whilst we find the contrary for the 
first group. The discrepancy in the first group may 
be attributed to wind clumping, but this would not 
explain the discrepancy for the objects below the BSJ 
temperature. An investigation into the effect of wind 
clumping on the predictions, for objects both above and 
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below the BSJ temperature, would shed light on this issue. 

Our radio study of 30 supergiants of spectral types 
08-B3 has revealed that for stars around the BSJ, the 
wind efficiency appear to be above the predicted general 
declining function towards later spectral types. The tem- 
perature of the potential local maximum also agrees with 
the predicted temperature of the BSJ. Possibly, a similar 
behaviour is present for the Ha results from SMC stars 
by Trundle k Lennon (2005). 

From the analysis of the WLR obtained here, we found 
the corresponding empirical CAK force multiplier to lie in 
the range Ues = 0.78 — 1. For these large values of a^s, 
the concept of the WLR may loose its meaning. Studies of 
the wind efficiency as a function of effective temperature 
should be considered a powerful alternative to the more 
often used WLR. 
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